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1000 



Select an initial atom (step 1010). 

+ . 

Assign a bond order to an adjacent bond (step 1015) 

Assign, to a next adjacent atom, an ESVD that is consistent with the 
previous bond (step 1020) 

Select the next adjacent bond's order (step 1030) 

Repeat the same approach for additional atoms (step 1040) 

y_ 

If an unacceptable state is encountered 



Reject the path values (step 1050) 

I 

Backtrack to ttie last point where a selection was made and proceed 
forward from there with a different selection (step 1060) 



I 



If a solution is found that cannot be improved upon, terminate (step 1070) 


u 




If a solution is found that is not optimal and there are other choices, 
explore the other choices (step 1080) 





FIG. 22 
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^2000 



Identify discrete (i.e., not adjacent) delocalized systems present in the 
structure (step 2010). 



Select any delocalized bond that has not already been treated (step 
2020) 



Add all adjacent delocalized bonds ^step 2030) 








Continue to add delocalized bonds adjacent to those just found, 
until no new bonds are added (step 2040) 



For each of the identified delocalized systems (step 2050) 



Analyze DS for characteristics (step 2060) 



Note the total charge and radical count required of the pi 
system (step 2070) 




^ 


r 




Calculate the internal coordination number ("ICN") of each 
atom (step 2080) 




r 


Identify bonds in DS that must be fixed, as a consequence of 
having an adjacent fixed multiple bond (step 2090) 



Develop a path through DS (step 2100) 
Construct the strategy list (step 2110) 

r 

Identify ring systems (step 2120) 
Set the best solution ("BS") to "undefined" (step 2130) 




FIG. 23A 
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9 

Execute until an identification of a solution deemed perfect, an 
exhaustion of possibilities, or an exceeding of an allocated amount 
of time (step 2140) 



Skip any mode that is incompatible with the control flags 
(step 2150) 

u 

Initialize the state (step 2160) 

I 

Queue the state ((step 2170) 



J 

While there are queued states, dequeue the most promising 
("S") and pass S to a dispatcher procedure (step 2180) 

y 

If the Dispatcher returns "Perfect Solution", skip to step 2210 



If BS is undefined after the mode processing is complete, return an 
indication of failure and terminate entire procedure (step 2200) 




r 




Apply BS to the given structure (step 2210) 





FIG. 23B 
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3000 







If the time elapsed has exceeded the amount of time allocated, return 
"Keep Going" (step 3010) 



Increment the strategy step index for S (step 3020) 




u 


Execute the procedure ("action") associated with the instant strategy step 
(step 3030) 


v' 


If the respective procedure returns "Perfect Solution", return "Perfect 
Solution" (step 3040). 


V 




If the respective procedure does not return "Perfect Solution" 






Decrement the strategy step index for S (step 3050 1 








Return "Keep Trying" (step 3060) 



FIG. 24 
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4000 



Use procedure 9000 to form a bitmask BM a that describes a (step 4010) 








Determine a set ("C") of ESVDs of a's element type that are compatible with 
a (step 4020) 



If no ESVDs qualify for inclusion in C, return "Failure" (step 4030) 



f 

Assign a penalty ("P") to each member (e) of C (step 4040) 



Set the penalty (P) to zero (step 4060) 

i 

If the ESVD represents a radical, and the state's residual radical flag 
is clear, increment P by three (step 4070) 



y 

If the ESVD is charged, and the residual charge of the state is non- 
zero and of opposite sign, add 50 to P (step 4080) 

J 

If the ESVD is charged but the state's residual charge is zero, 
increment P by 2 (step 4090) 

j 

If e is charged and a is carbon, increment P (step 4100) 



f 

If e does not have an internally directed multiple bond, increment 
P (step 4110) 

t 

Sort the members of C in order of increasing penalty (step 4120) 

FIG. 25A 
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Queue inferior ESVDs for later consideration (step 4130) 



For each member (C.) of C except the first (i.e., the best), execute 
the following steps 4150 - 4180 (step 4140) 



Copy the current state is copied to a new state, q (step 4150) 
Set the ESVD of atom a in q to C. (step 4160) 

i 

Increment the strategy step counter for q (step 4170) 
Insert q into the priority queue (step 41 80) 



Apply the best ESVD to S (step 4190) 



Update the state variables having to do with atomic properties 
(step 4210) 



Subtract the ESVD's charge from the state's residual charge 
(step 4220) 

i ~ 

If the ESVD has an unpaired electron, toggle the state's 
residual radical flag (step 4?m\ 



1 




Invoke the Dispatcher recursively (step 4240) 




u 




If the Dispatcher returns "Perfect Solution", return the same and terminate 


(step 4250) 







Restore S to its value before the best ESVD was applied (step 4260) 




FIG. 25B 
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1 




Return "Keep Tn/mg" (step 4270) 



FIG. 25C 
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5000 



Calculate the ASB of the preceding atom 


(step 5020) 




Calculate the difference BM diff = BM ESVD - 


BM, (step 5030) 


u 



Determine the possible bond orders (step 5040). 



Set AddSingle to true if and only if BM^has either the internal single 
bit or the two internal single bit set (step 5050) 



J 

Set AddDouble to true if and only if BM dijf has the internal double bit 
or the two internal double bit set (step 5060) 

f 

Set CooptExtern to true if a's ICN is three or greater, and a's valence 
allows a one or more units of bond order beyond the units already 
used up by a's fixed internal bonds (step 5070) 

* ~ 

If either AddSingle or CooptExtern is set to true 



Set b's order to single and any relevant state variables are updated 


(step 5100) 




v 




Invoke the Dispatcher recursively (step 5110) 




v 




If the Dispatcher returns "Perfect Solution", return the same (step 


5120) 





FIG. 26A 
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f 



If AddDouble is set to true 




Set b's order to double (step 5130) 








Invoke the Dispatcher recursively (step 5140) 




v 




If the Dispatcher returns "Perfect Solution", return the same (step 
5150) 


i 



Otherwise, return "Keep Trying" (step 5150) 



FIG. 26B 
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Calculate a's ASB using procedure 9000 



Correct BM a for tertiary atoms (step 6030) 



If BAi equals {0,1,2}, set BM a {0,2} (step 6040) 

* 

It the ESVD calls for U units of external bonding, but the atom does not 
express U units of external bonding, and the control flags permit implicit 
hydrogens, include bit 4 in BM a (step 6050) 









Recursively invoke the Dispatcher (step 6060) 












Return whatever value the Dispatcher returns (step 6070) 










i 



Otherwise, return "Keep Trying" (step 6080) 



r 



6000 



(step 6020) 



FIG. 27 
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7000 



If S has a residual charge or radical, return "Keep Trying" (step 7020) 








If mesomers are being enumerated, offer the dekekulized structure back 
to source that invoked dekekulization (step 7030) 


* 




Calculate S's rating (R) (step 7040) 


I 


If R is zero and mesomers are not being enumerated, return "Perfect 
Solution" (step 7050) 








If the best solution is undefined, or if R is better than the rating of BS, set 
BS to S (step 7060) 


1 




Return "Keep Trying" (step 7070) 



FIG. 28 
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8000 



Initialize ESB to zero (step 8010) 

I . 

Set bits of ESB according to the ESVD's charge (step 8020) 
Set bits according to the ESVD's radical character (step 8030) 



Set bits according to the ESVD's pair of internal bond orders (step 8040) 






r 




Set bits according to whether the ESVD has one or more units of external 
bonding (step 8050) 



Return ESB (step 8060) 



FIG. 29 
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9000 



Initialize ASB to zero (step 9010) 



J 

Set bits of ASB according to the number of internal bonds adjacent to a that 
are fixed, with bond order single (step 9020) 



J 

Set bits according to the number of internal bonds adjacent to a that are 
fixed, with bond order double (step 9030) 



t 

If a has more than two adjacent bonds, set bit {4} (step 9040) 

* 

If a has exactly two adjacent bonds, and if the control flags dictate that no 
implicit hydrogens exist, and a's element type is not carbon, set bit {5} 
(step 9050) 

y 

If the residual radical value of S is "no-radical", and the control flags do 
not permit unnecessary creation of radicals, set bit {14} (step 9060) 

y 

If the residual charge of S is zero and the control flags do not permit 
unnecessary creation of charged atoms; or if a's element type is carbon 
and the control flags require charges to be situated on heteroatoms, set 
bits {11,12} (step 9070) 



v 



Otherwise 






If the residual charge of S 


is positive, set bit {11} (step 9075) 












If the residual charge of S 


is negative, set bit {12} (step 9080) 













Return ASB (step 9090) 



FIG. 30 
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Initialize R to zero (step 10010) 

I 

For each ring system ("RB") in DS, execute steps 10030 - 10250 (step 10020) 



Set periphery (P) to "empty" (step 10030). 

I 

For each ring r belonging to the RB, execute steps 10050 - 10110 
(step 10040) 



Set P to the exclusive OR of itself with the bonds in r (step 
10050) 



Assess the "one-ring" penalty as follows (steps 10070-10080) 
(step 10060) 




If e is a multiple of four ("4n") and the control flags 
prescribe the penalization of anti-aromatic solutions, 
subtract two from R (step 10070) 












If e is odd, subtract one from R (step 10080) 







For each ring r2 belonging to RB wherein r2 > r, execute 
steps 10100 - 10110 (step 10090) 



If r and rl have one or more rings in common 



Set RC to be the compound ring of r and rl 
(step 10100) 



J 

Adjust R by the one-ring penalty amount, 
where rl is substituted for r (step 10110) 



^10000 



FIG. 31A 
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If there are three or more rings in RB 




Take r3 to be the compound ring represented by P 
(step 10120) 








Adjust R by the one-ring penalty amount, where r3 is 
substituted for r (step 10130) 



Apply a significant penalty for gratuitous charges and radicals 
(steps 10160 - 10180) (step 10140) 



Take NSC to be the number of superfluous charges (step 


10160) 






r 


Take NSR to be the number of superfluous radicals (step 


10170) 






f 


Reduce the score by ten times (NSC + NSR) (step 10180) 



Apply a small penalty tor locating a charge on a carbon rather than 
an available heteroatom (steps 10200 - 10220) (step 10190) 



Take NCC to be the number of charged carbon atoms in S 
(step 10200) 



I 



Take NNH to be the number of neutral heteroatoms in S 
(step 10210) 




V 




If NCC exceeds NNH, reduce R by (NCC - NNH) (step 10220) 



Penalize pairs of adjacent atoms having identical, either empty or 
full, orbitals (step 10230) 

FIG. 31B 
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9 

Return R (step 10250) 



FIG. 31C 



